The effects of violating detailed balance on critical dynamics 
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We present an overview of the effects of detailed-balance violating perturbations on the universal 
static and dynamic scaling behavior near a critical point. It is demonstrated that the standard crit- 
ical dynamics universality classes are generally quite robust: In systems with non-conserved order 
parameter, detailed balance is effectively restored at criticality. This also holds for models with con- 
served order parameter, and isotropic non-equilibrium perturbations. Genuinely novel features are 
found only for models with conserved order parameter and spatially anisotropic noise correlations. 
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One of the major goals in theoretical non-equilibrium 
physics has been the identification and classification of 
universality classes for the long- wavelength, long-time 
scaling behavior both near continuous dynamical phase 
transitions, and for systems displaying generic scale Ln- 
variance. Indeed, through investigations of certain spe- 
cific models, a number of prototypical non-equilibrium 
universality classes have been identified. Prominent ex- 
amples are driven diffusive systems 0j, models of driven 
interfaces and growing surfaces Q , depinning transitions 
[3|, and phase transitions from active to absorbing states 
\i |, e.g., in diffusion- limited chemical reactions. 

A complementary approach is to study the influence 
of non-equilibrium perturbations on the known univer- 
sality classes for equilibrium dynamical critical phenom- 
ena p. Equilibrium critical dynamics is concerned with 
the relaxational and reversible kinetics near a thermody- 
namic critical point at temperature T c , as generically de- 
scribed by the Landau- Ginzburg- Wilson (LGW) model 
for an n-component order parameter vector field S in 
d space dimensions p|. In addition to the two inde- 
pendent static critical exponents, e.g. the correlation 
length exponent v defined via £ cx \t\~ u (t = T — T c ) 
and Fisher's exponent r\ for the algebraic decay of the 
two-point correlation function at criticality (T = T c ), 
C(x — x') cx | a; — x '\~( d ~ 2+ri \ the order parameter relax- 
ation is governed by a dynamic exponent z that describes 
critical slowing down: The characteristic time scale di- 
verges as < c h cx \t\~ zv upon approaching the transition. 
This allows for time scale separation and thus a formula- 
tion of critical dynamics in terms of non-linear Langevin 
equations: The relevant 'slow' modes consist of the order 
parameter and all conserved quantities to which it is stat- 
ically or dynamically coupled. All remaining 'fast' de- 
grees of freedom are captured through an effective Gaus- 
sian white noise. Different values for z ensue depending 
on whether the order parameter is a conserved quan- 
tity or not, and on the additional conserved quantities 
present. The diffusive relaxation of the latter near criti- 
cality can either be characterized by the same exponent 
z as for the order parameter ('strong' dynamic scaling), 



or be given by different power laws ('weak' scaling) [pi. 

In order to ensure relaxation towards thermal equilib- 
rium at long times, as given by a Gibbs distribution, one 
has to carefully implement detailed-balance conditions. 
In the language of non-linear Langevin equations, these 
are (i) the Einstein relation between the relaxation con- 
stants and the noise strengths, and (ii) the condition that 
the probability current associated with reversible kinetics 
be divergence- free @. Naturally, the following question 
arises: What happens if the equilibrium conditions are vi- 
olated? In this letter, we shall explore two generic types 
of detailed-balance violations, namely (a) coupling the 
order parameter and additional conserved quantities to 
heat baths with different temperatures, and (b) allowing 
for spatially anisotropic noise correlations for conserved 
variables. To determine their universal features, we map 
the Langevin dynamics for most of the models listed in 
Ref. H to a dynamic field theory ||], and employ stan- 
dard renormalization group (RG) methods . 

In the theory of static critical phenomena, an analo- 
gous issue concerns the effect of terms that break the 
original order parameter symmetry, e.g., the influence of 
cubic anisotropies on the isotropic n-component Heisen- 
berg model. There, the rotational symmetry is restored 
at criticality, provided n < n c w 4 ||. Since de- 
tailed balance originates from time-reversal symmetry, 
we might anticipate that it could, under certain condi- 
tions, effectively become reinstated in non-equilibrium 
critical dynamics, whereupon the asymptotic scaling laws 
are those of the corresponding equilibrium model. Yet 
even then, the question arises whether there exist any 
non-equilibrium dynamical RG fixed points that could 
strongly influence crossover regimes. The other possi- 
ble scenario is of course that violating detailed balance 
constitutes a relevant perturbation, rendering the equi- 
librium RG fixed point unstable, and driving the system 
towards a genuine non-equilibrium stationary state. 

The simplest case represents purely relaxational dy- 
namics, with either non-conserved or conserved, and 
therefore diffusively relaxing, order parameter S |1C| ] 
(models A and B according to the classification in Ref. [p|, 
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respectively). With the effective LGW Hamiltonian 



d t S(x, t) = X ± V 2 ± 5H cS [S]/5S{x, t) + C(x, t) 



(0) 



H[S] 



d d x 



.S(x) 2 + i[VS(x)] 2 + ^S(xY 



(1) 



the corresponding Langevin equations read 



dtS(x, t) = -X (iV) a 5H[S]/SS(x, t) + ((x, t) (2) 

where a = (2) for models A (B). The noise has zero 
mean, (£(x, t)) = 0, and its correlations are taken to be 

(C(x,t)C {x',t')) = 2\(W) a 5 al3 6(x-x')5(t-t') . (3) 

The equilibrium Einstein relation requires A = k^T A, 
which ensures that the associated probability distri- 
bution for a configuration S approaches P e q[S] oc 
cxp(— i/[S]/fceT) as t — > oo. Setting A ^ A violates 
detailed balance; yet this simply represents a tempera- 
ture shift here. Indeed, upon rescaling S — * (A/A) 1 / 2 S 
and u — ► u = (A/A) 1 / 2 ^, detailed balance is formally re- 
stored, albeit with a different non-linear static coupling 
jyj . At criticality, the RG flow will take the latter to the 
Heisenberg fixed point, independent of its initial value, 
and thus the equilibrium critical exponents are recov- 
ered, namely to lowest non-trivial orders in the devia- 
tion e = 4 — d from the upper critical dimension d c = 4: 
v- 1 = 2-(n+2)e/(n+8), n = (n+2)e 2 /2(n+8) 2 @, and 
z = 2 + en with c = 61n| — l+0(e) for model A, whereas 
z = 4 — n (exactly) for model B Jul . We already note 
that the same rescaling with g 2 — »lp — Xg 2 /X maps the 
isotropic 0(3)-symmetric non-equilibrium model J 

d t S = -g SH[S]/SS xS + AV 2 SH[S]/5S + ( , (4) 
(C(x,t)C (x',t')) = -2XW 2 S a0 S(x-x')S(t-t') (5) 

onto its equilibrium version, with z = (d + 2 — f])/2 [ jl2| . 

It is indeed a well-established fact that model-A or 
Glauber dynamics for the kinetic Ising model (n = 1) 
is quite robust against non-equilibrium perturbations 
[13,[l4|], even when these break the up/down symmetry 
[ 1 5 1 . Novel features only arise when Kawasaki dynam- 
ics is introduced, whereupon the order parameter be- 
comes conserved (model B), and in addition the noise 
strength is rendered anisotropic: AV 2 — * A|| V 2 + A^V^ 
in Eq. (J3J) , corresponding to coupling the longitudi- 
nal and transverse sectors (of dimensions d\\/±) to heat 
baths with different temperatures Tiuj_. For this two- 
temperature or randomly driven model B, and for Tu > 
Tj_, only the transverse sector softens at T c . As a con- 
sequence, the system develops marked anisotropics akin 
to driven lattice gases (or equilibrium uniaxial dipolar 
magnets, Lifshitz points, and ferroelastic materials), de- 
scribed by the generalized scaling law C(q\\ , qx, v, r) = 
qj2-*+*C±(qxt,q\\/q? A ,u/q ± ) §. As A ^ 0, equilib- 
rium isotropic scaling is violated. However, upon omit- 
ting irrelevant terms, the Langevin equation at criticality 
may be written in the form 



with an effective Hamiltonian that contains long-range 
interactions fTl [here, / . . . = / d d q/(2n) d . . .] 



cq 2 + q 2 x(rx + q 2 x) 
2^ 



\S(q)\ 2 + ld d x^& 



(J) 



These strongly affect the nature of the ordered phase, 
and reduce the upper (to d c = 4 — du) as well as lower 
critical dimension p6jl7 |. To one- loop order (only), the e 
expansion for the static critical exponents formally yields 
the same results as for the equilibrium Heisenberg model, 
albeit with e = 4 — d— d\\. The two-loop results for d\\ — 1 
can be found in Ref. p8| . Moreover, the exact scaling 
relations A = 1 — r)/2 and z = 4 — 77 hold. 

We have extended these considerations to models C 
and D ]l9|], which take the static coupling 



AH[S,p} = J d d a 



l -p{x) 2 + g -p{x)S{x) 2 



(8) 



of the non-conserved/conserved order parameter to the 
conserved energy density p into account pcfl . The dy- 
namics is then defined through Eq. (||) with H[S,p] — 
H[S] + AH[S,p], setting a = 0(2) for model C (D), and 

d t p(x, t) = DV 2 6H[S, P ]/5p{x, t) + V (x, t) , (9) 
(r)(x, t) r](x',t')) = -2D V 2 5{x - x')S(t - t') . (10) 

In addition to u, there appears now a three-point cou- 
pling / = Xg 2 /X, also marginal in d c — 4, as well as 
the dimensionless ratio of time scales w = X/D, and 
the parameter = DX/DX, which can be interpreted 
as the temperature ratio T p jTs of the heat baths cou- 
pled to the energy density and order parameter, respec- 
tively. For model C, we find the associated RG one-loop 
beta function e = -26(1 - 9) / w/{l + w) 2 . This re- 
sult already establishes the stability of the equilibrium 
fixed point 0* q = 1. In fact, we have not found any 
other, genuinely non-equilibrium RG fixed points cor- 
responding to O = or = 00 at all. This remains 
true even when we allow for spatially anisotropic noise 
for the conserved field p, as the RG beta function for 
the additional new variable a = 0||/0j_ becomes to 
one-loop order: f3 a = —a(l — <t)/0_l/3(1 + w) 3 , with 
the isotropic fixed point a* s = 1 being stable. Thus, 
there are merely the three (one-loop) equilibrium scal- 
ing regimes pC( |: (a) Strong-scaling, n = 1: w* = 1, 
zg = z p = 2 + a/v; (b) weak-scaling, 2 < n < 4: w* = 00, 
z s = 2(1 + a/nv) < z p = 2 + a/v; (c) 'model A', n > 4: 
w* = 00, zs = 2 + ci], z p = 2; here a — 2 — dv denotes the 
specific-heat critical exponent. In contrast to the above 
results, a linear coupling of a conserved mode to the or- 
der parameter induces effective long-range interactions, 
leading to novel genuine non-equilibrium behavior |2lJ|. 
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For model D, the characteristic time scale for the or- 
der parameter is always much larger than that of the 
diffusive field, and therefore w* = 0. By integrating out 
the scalar fields p from the dynamic functional, one can 
then readily show that the non-equilibrium parameter 
disappears entirely from the theory in the limit w — > 
p9[ , which completely reduces the model to the equilib- 
rium theory, with zs = 4 — rj, whereas z p = 2 + ajv for 
a > and z p — 2 for a < ('model B'). Upon allowing 
for spatially anisotropic order parameter noise correla- 
tions, we arrive at a two-temperature model D. As for the 
corresponding pure order parameter relaxation dynam- 
ics, the non-equilibrium perturbations then induce strong 
anisotropies, whereupon the effective Hamiltonian (0) en- 
ters the model-D equations of motion, with only trans- 
verse Laplacians, and the associated downwards shift of 
d c . The model-D scaling relations for the anisotropy and 
dynamic exponents still apply, as well as the two different 
regimes for a > and a < 0, albeit with the static ex- 
ponents v, 77, a of the long-range anisotropic theory jl9) . 

Purely relaxational dynamics provides, however, an in- 
sufficient description for many real systems. Often, fully 
reversible, non-dissipative terms originating in the mi- 
croscopic dynamics need to be taken into account; in 
isotropic ferromagnets (model J), for example, there is 
also the spin precession in the local magnetic field, see 
Eq. (Q). Such reversible contributions in the Langevin 
equations for critical dynamics may also involve mode- 
couplings to other conserved and therefore slow variables: 
In planar ferromagnets, the (non-conserved) order pa- 
rameter is confined to the xy plane, say, but the non- 
vanishing commutators of the spin components yield a 
coupling to the diffusive fluctuations of the conserved z 
component of the magnetization ((S z ) = 0). In isotropic 
antiferromagnets, the order parameter is represented by 
the three-component staggered magnetization (not con- 
served), dynamically coupled to the conserved magne- 
tization. The corresponding Langevin equations define 
models E and G, respectively [p2| ; their generalization to 
n order parameter components is termed the SSS model, 
originally introduced in the context of structural phase 
transitions [^3],|4| . Lastly, a consistent description of the 
critical dynamics near the liquid-gas transition, or equiv- 
alently, of the phase separation in binary liquids, involves 
not only the conserved scalar order parameter density 
S (a linear combination of the mass and energy densi- 
ties) , but in addition the independent and also conserved 
transverse corresponding current density jj_ (model H) 
P^,|3. In the following, we describe the effects of non- 
equilibrium perturbations on these dynamic universality 
classes with reversible mode-couplings . 

The results of our investigations of isotropic and 
anisotropic non- equilibrium versions of the SSS model 
were already reported in Ref. pJ] ; for completeness we 
review our essential findings here. As in model C, we 
may choose different effective temperatures for the n- 



component order parameter S and the n(n — l)/2 non- 
critical, conserved generators M of the rotation group 
0(n). The coupled Langevin equations of motion read 

d t S a = gJ2M af3 SP -\5H[S]/5S a + C , (11) 

{( a (x,t)(P(x',t')) = 2X5 at3 S(x-x')S(t~t r ) , (12) 
d t M af3 = -g (S a V 2 SP - S p \7 2 S a ) + DX7 2 M af3 + rj af3 , 

(13) 

(n af) (x, t) rf s {x', t')) = -2DV 2 8(x - x')5{t - t') 

x {5 af3 5~i s - S^S^ 1 ) . (14) 

As opposed to the purely relaxational models, there exist 
genuine non-equilibrium fixed points here, correspond- 
ing to temperature ratios = Tm/Ts — and 00. In 
the former case, there is no feedback to the order pa- 
rameter, and the critical dynamics is essentially model- 
A like (zs = %m = 2), with anomalous noise corre- 
lations oc q d ~ 2 for the generators M. For 0^ = 00, 
the order parameter dynamics does not affect the con- 
served fields, thus zm — 2, but is itself strongly influ- 
enced by their fluctuations (w = X/D = 00). We find 
zs = d/2, as in equilibrium, but a modified static expo- 



nent v~ = 2 — 2(n + 2)e/(n + 



111. However, both 



these fixed points are unstable, and the asymptotic crit- 
ical properties are governed by the equilibrium strong- 
scaling fixed point with zs = Zm — d/2 p3| , p4| . Allow- 
ing for spatially anisotropic noise for the generators M 
provides the new parameter a = 0||/O_l, but does not 
change the overall picture: The isotropic fixed point o~* s 
remains stable, whence the equilibrium critical behavior 
is eventually recovered . The anisotropic fixed points 
are essentially determined by a combination of the = 
and = 00 characteristics in the different spatial sectors. 

The intriguing question now arises how non- 
equilibrium perturbations affect the critical dynamics of 
a conserved order parameter field, when reversible mode- 
couplings are present. We have therefore investigated the 
non- equilibrium model H, as defined by p6[ 



d t S = ~g{VS) • j + AV 2 t - V 



^ 2 )* + c 



(C(x, t) ((x',t')) = -2A V 2 S(x - x')S(t - t 1 ) , 
d t h = T [g(VS) (t - V 2 ) S + D\7 2 i + v ] , 



(15) 

(16) 
(17) 



(77(2, t) ■ t](x', t')) = -2dD V 2 S(x - x')S(t - t') , (18) 

where T(q) = Sij — qiqj/q 2 denotes the transverse pro- 
jector in momentum space, and we have omitted irrele- 
vant (in the RG sense) terms. As in model D, asymp- 
totically w — X/D — * 0; as a consequence, only cme 
additional non-equilibrium fixed point 0^ = is al- 
lowed. In analogy with the corresponding SSS model 
fixed point, the order parameter dynamics is model-B 
like, zs = 4, while Zj = 2, with anomalous noise correla- 
tions again oc q d ~ 2 in Fourier space Eg]. Yet from (}q = 



3 



-(1 - 9)§/(l + 1/16 8), where / = \g 2 /\ 2 D, we infer 
that once more the equilibrium weak-scaling fixed point 
Q* E = 1 is stable, with z s = 4 - 18e/19 and zj = 2 - e/19 
(Galilean invariance fixes zs + Zj = d + 2) pq , ^4] . How- 
ever, the generalization to spatially anisotropic noise, 
A V 2 -» Ay V| + A ± Vi and D V 2 -> D\\ Vf + D±_ V\ 
has a drastic effect. As in the corresponding models B 
and D, in this two-temperature model H, the characteris- 
tic anisotropic DDS singularities emerge. In combination 
with the mode-coupling terms, this prevents the system 
from approaching an effective equilibrium model: The 
divergence-free condition for the reversible probability 
current cannot be satisfied. In fact, to one-loop order 
at least we even find a run-away flow for the RG cou- 
plings, and are thus unable to determine the long-time 
scaling behavior f26fl . Remarkably, the similarly con- 
structed anisotropic, two-temperature model J is plagued 
by the same pathology as is a recently studied uni- 
formly rather than randomly driven model J pjj] . In that 
instance, computer simulations revealed that the system 
displays spatio-temporal chaos at long times; perhaps 
the absence of an RG fixed point in the two-temperature 
models J and H may indicate chaotic behavior as well. 

We have investigated the effect of detailed-balance vi- 
olations on critical dynamics. Generally, models with 
non-conserved order parameter are quite robust against 
non-equilibrium perturbations. The relaxational models 
A and C only have an equilibrium fixed point. For the 
SSS model, comprising models E and G, genuine non- 
equilibrium fixed points do exist, corresponding to uni- 
directional couplings between order parameter and con- 
served fields, but are unstable. Thus, at criticality, the 
standard critical behavior is eventually recovered. This 
remains true even when the conserved noise is rendered 
spatially anisotropic. Essentially the same statements 
apply for models B, D, J and H with conserved order 
parameter, provided detailed-balance violations are in- 
troduced isotropically. With spatially anisotropic order 
parameter noise correlations, however, we find (to one- 
loop order) no RG fixed points for models J and H with 
reversible mode-coupling terms. In contrast, the two- 
temperature relaxational models B and D are asymptot- 
ically described by an effective equilibrium model, with 
characteristic anisotropic, long-range correlations. 
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